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A  time -domain  boundary  integral  equation  method  has  been  developed  to 
calculate  elastodynamlc  fields  generated  by  the  incidence  of  stress  (or 


displacement)  pulses  on  single  cracks  and  systems  of  two  collinea^cracks . 
The  system  of  boundary  integral  equations  has  been  cast  in  a  form  which  is 
amenable  to  solution  by  the  boundary  element  method  in  conjunction  with  a 
time-stepping  technique.  Particular  attention  has  been  devoted  to  dynamic 
overshoots  of  the  stress  Intensity  factors.  Elastodynamlc  stress  intensity 
factors  for  pulse  incidence  on  a  single  crack  have  been  computed  as 
functions  of  time,  and  they  have  been  compared  with  results  of  other 
authors .  For  collinear  macrocrack-microcrack  configurations  the  stress 
Intensity  factors  at  both  tips  of  the  macrocrack  have  been  computed  as 
functions  of  time  for  various  values  of  the  crack  spacing  and  the  relative 
size  of  the  microcrack. _ 
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The  interaction  of  a  stress  pulse  with  a  crack  may  give  rise  to  dynamic 
overshoots  of  the  stress  intensity  factors,  i.e.,  larger  stress  intensity 
factors  than  would  be  obtained  for  quasi-static  application  of  the  external 
loads.  For  step-stress  pulses  the  dynamic  overshoot  phenomenon  has  been 
noted  by  Achenbach  [1],  Thau  and  Lu  [2]  and  Sih  [3].  For  the  two 
dimensional  configuration,  Thau  and  Lu  [2]  used  integral  transform 
techniques  in  conjunction  with  a  generalized  Wiener-Hopf  method  to  obtain 
exact  but  short-time  dynamic  stress  intensity  factors.  Sih  et  al.[3] 
applied  integral  transform  methods  together  with  a  numerical  solution  of 
dual  integral  equations  and  a  numerical  inversion  of  Laplace  transforms  to 
calculate  dynamic  stress- intensity  factors.  The  problem  has  also  been 
treated  by  finite -element  [4]  and  finite  difference  methods  [5].  An 
extension  of  the  approach  of  Ref. [3]  has  been  presented  by  Itou  [6],  who 
analyzed  the  fields  for  transient  wave  interaction  with  two  coplanar  cracks 
of  equal  length. 

In  this  paper  we  investigate  pulse -generated  crack- tip  fields  by  the 
use  of  a  time-domain  boundary  Integral  equation  (BIE)  method.  For  wave 
interactions  with  volume  scatterers,  this  technique  has  been  successfully 
applied  to  two-dimensional  elastodynamic  problems  by  several  authors,  see, 
e.g.,  Niwa  et  al.[7],  Manolis  [8],  Mansur  [9],  Antes  [10], [11],  and  Estorff 
[12] .  The  usual  displacement  BIE  formulation  for  scattering  of  elastic 
waves  by  volume  scatterers  degenerates,  however,  when  the  scatterer  is 
reduced  to  a  flat  crack.  A  remedy  for  this  difficulty  is  the  use  of 
"traction"  BIE's,  which  results  in  a  system  of  singular  integral  equations 
for  the  unknown  crack  opening  displacements.  Unfortunately,  such  BIE’s  are 
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highly  singular  and  they  cannot  he  solved  directly  by  numerical  methods.  In 
the  present  paper  this  difficulty  is  overcome  by  reducing  the  higher  order 
singularities  to  integrable  singularities  which  can  be  incegraced 
analytically  or  numerically.  The  simplified  time -dependent  BIE's  are  solved 
by  the  boundary  element  method  in  conjunction  with  a  time -stepping 
technique.  An  alternative  method  has  been  proposed  by  Nishimura  et  al.[13] 
who  used  a  double  layer- regularization  procedure.  The  numerical  approach  of 
the  present  paper  has  been  applied  to  obtain  time  histories  of  elastodynamic 
stress  intensity  factors  for  single  cracks  as  well  as  for  configurations  of 
a  macrocrack  and  a  collinear  microcrack.  Parametrical  studies  show  the 
influence  of  the  size  and  location  of  the  microcrack  on  the  effective  stress 
intensity  factors  of  the  solitary  crack. 

The  results  for  collinear  cracks  approximate  the  fields  for  macrocrack- 
microcrack  configurations  that  are  often  observed  in  brittle  materials  such 
as  ceramics,  rocks  and  concretes.  In  such  solids  the  high  level  of  stress 
and  deformation  in  the  vicinity  of  a  crack  tip  gives  rise  to  microcracking 
and/or  the  formation  of  microvoids  in  a  confined  zone  surrounding  the  macro - 
crack  tip  (see  Bradt  et  al.[14],  Carpinteri  et  al.[15],  Hoagland  et 
al.[16]).  The  existence  of  neighboring  micro-cracks  may  significantly  alter 
the  stress  intensity  at  the  main  crack  tip,  as  shown  by  Kachanov  et  al.[17], 
Rubinstein  [18]  and  Rose  [19]  for  the  static  case.  Depending  on  the  size 
and  location  of  microcracks  or  microvoids,  their  presence  can  either 
increase  the  stress  intensity  factors  (stress  amplification)  or  decrease  it 
(stress  shielding  or  toughening).  Knowledge  of  the  dependence  of  the  stress 
intensity  factors  on  the  microdefects  will  assist  in  predicting  macrocrack 
propagation.  For  static  loading  several  studies  can  be  found  in  the 


literature,  see  Chen  [20],  Chudnovsky  et  al . [21] , [22] ,  Kachanov  [23], 
Kachanov  et  al.[17].  Rose  [19],  Rubinstein  [18] , [24]  and  Yokobori  et 
al. [25] . 

Recently,  Zhang  and  Achenbach  [26], [27]  have  investigated  the  effects 
of  microvoids  and  microcracks  for  time-harmonic  wave  loading.  It  was  found 
that  dynamic  effects  may  give  rise  to  considerable  amplifications  of  the 
stress  intensity  factors. 


A  two-dimensional  configuration  of  a  homogeneous,  isotropic,  linearly 
elastic  body  containing  a  macro-crack  and  a  collinear  neighboring  micro¬ 
crack  is  shown  in  Fig.  1.  An  incident  displacement  pulse  interacts  with  the 
two  cracks,  and  generates  a  scattered  displacement  field.  The  propagation 
direction  of  the  incident  pulse  is  in  the  Xj^X2-plane,  and  the  analysis  of 

this  paper  is  two-dimensional  and  for  a  state  of  plane  strain.  In  terms  of 

in  sc 

the  incident  field,  u  (x,t),  and  the  scattered  field,  u  (x,t),  the  total 

a  —  a  — 

displacement  field  may  be  written  as 


u  (x,t)  -  u^"(x,t)  +  u®‘^(x,t)  ,  a  -  1,2. 
o  —  a  —  a  — 


(2.1) 


Similarly  we  have  for  the  stress  components 


(2.2) 


Since  the  faces  of  the  cracks  are  free  of  tractions,  the  following 
conditions  follow  for  the  scattered  field 


(x,t) 


for  X  «  +  Fj  , 


(2.3) 


where  and  F2  define  the  faces  of  the  macro-crack  and  the  micro-crack, 
respectively.  The  initial  conditions  are 


u^^(x,t)  -  u®‘^(x,t)  -  0  for  t  <  0. 

Qt  —  O  — 


(2.4) 


Here  time  t  starts  when  the  incident  wave  first  reaches  the  crack. 

The  integral  representation  for  the  components  of  the  scattered 
displacement  field  may  be  written  as 


sc 

u 
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(Xp.t) 


(2.5) 


Q 

Here  Xp  is  the  position  vector  of  the  observation  point,  0^2^  stress 

Green's  function  for  the  uncracked  plane  (see  Appendix  A),  and  Au^  is  the 
displacement  jump  across  the  crack  faces  (the  crack- opening  displacement) , 


Au^  (x^,0,r)  -  u^(Xj^,0  ,r)  -  u^(x^,0  ,r) 


(2.6) 


As  shown  in  the  next  section,  the  solution  to  the  problem  formulated  in  this 
section  can  be  reduced  to  the  solution  of  a  set  of  boundary  integral 
equations . 


wr*  vv  ■.-*  ^.-t  (.-y -.-v  iv  ^ 


3 .  Derivation  of  Discretized  BIE's 
In  the  conventional  procedure  BIE's  can  be  derived  from  Eq.(2.5)  by 
taking  the  limit  Xp  -*  F  -  +  r2.  However,  such  "displacement"  BIE's 

degenerate  for  cracks ,  and  hence  they  are  not  a  valid  basis  for  numerical 
modeling.  This  difficulty  is  overcome  by  the  use  of  "traction"  BIE's. 
These  are  obtained  by  substituting  Eq.(2.5)  into  Hooke's  law 


-  XS  „\i  + 

7,7 


a,P 


p,a 


(3.1) 


which  yields  a  representation  formula  for  the  stress  components  0^2 
observation  point  Xp  as: 


o  F 


(3.2) 


where  n^  are  the  components  of  the  normal  vector  to  F,  and 
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(3.3) 


Then,  BIE's  can  be  obtained  from  Eq.(3.2)  by  taking  Xp  -►  F.  Unfortunately, 
such  "traction"  BIE's  are  hyper  singular  when  the  observation  point  Xp  and 

Q 

the  source  point  x  coincide,  since  in  this  case  the  kernel  function  o 

Eq.(3.2)  behaves  as 


if  r  -►  0 


(3.4) 


and 


q(t-r)2-rVc2]® 


,  if  (t-r)2  -►  r*/c2 


(3.5) 


in  which  r  -  x-x  ,  and  c  is  either  c.  or  c_,  where 

•  «.  .-p  ‘  Irf  1 


c2  -  (A+2/i)/p  ,  c2  -  n/p. 


(3.6) 


A  detailed  discussion  of  these  singularities  can  be  found  in  a  paper  by 
Nishimura  et  al.  [13]. 

To  reduce  these  higher  order  singularities,  Nishimura  et  al.  [13] 
proposed  a  double  layer -regularization  procedure.  In  the  present  paper  we 
apply  another  regularization  method  which  has  been  developed  by  the  authors 
in  a  recent  paper  [26]  for  scattering  of  incident  time -harmonic  elastic 
waves  by  multiple  cracks. 

Following  the  procedure  of  Zhang  and  Achenbach  [26],  we  first  divide  F 
into  J  elements.  Then  Eq.(3.2)  can  be  written  in  the  following  discretized 


J  t  s 

ff®2^Xp,t)  -  -  S  /  J  K^25e(2p-*^:5>’’)^’^5(2S-’’)^jdsdT, 


j-1  o  s. 


2p  ^  r 


(3.7) 


where  Sj  and  tihe  endpoints  of  the  j-th  element.  Because  the  main 
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structure  of  the  three  terms  in  (see  Eq.(3.3))  is  similar,  we  will 

consider  only  the  following  integral 


I,  -  J  ■'  <7^ 

j  i  Sey.y  6  « 

j 


(3.8) 


By  adding  and  subtracting  the  same  terms,  Eq.(3.8)  can  be  rewritten  as 


I,  -  J  (<7?  Au^n  -  at  Au<.n  )ds 
j  J  S<7,7  5  «  fi«7.e  S  7 

j 


+  r  Au.n  ds 

i  Sey.e  S  y 


It  can  be  shown  that  the  first  integral  of  (3.9)  can  be  put  into  the 
following  form  (see  Appendix  B  of  Zhang  and  Achenbach  [26]) 


t  Of  AUf 

ye  eo7  d 


j+1  s 


r  j+l  G  ^ 

J  t  a  ,  Au,.  ,n  ds 


s . 
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where  e^^  denotes  the  two-dimensional  permutation  tensor.  The  second 


integral  of  Eq.(3.9)  can  be  rewritten  as 


r^j+l  G  .  r®j+l  . 

J  ■’  Or  Au^n  ds  -  p  J  d-  Au^n  ds 
J  5e7,e  5  7  J  §7  5  7 


where  the  equation  of  motion  for  the  Green' s  function 


has  been  used.  Here  is  the  Green's  function  for  the  displacement 


components  (see  Appendix  A).  Equation  (3.10)  and  Eq.  (3.11)  together  result 
in 


I.  -  e  a  .  AUf 
J  7<  <<^7  0 
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s . 
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^  r®j+l  ,G  ^ 

+  p\  'If  Au.n 
07  5  7 


s . 
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ds 


(3.13) 


Using  the  same  idea  described  above  for  the  second  and  the  third  terms  of 
Eq.(3.3)  we  obtain 
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T  ^ at  A'AfO  ds  -  H'^j-Au 

•1  apS  (St  aPS 
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r^j+l  HI  A  fi 
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(3.14) 


where 


“i^s 


G  G  G 

A5  .(  Oc  +  Pi^a  +  £  a-  q) 

ap  7£  Sty  Pe  Sta  at  StP 


(3.15) 


H2 

apS 


(3.16) 


Furthermore,  it  can  be  shown  that  the  following  relation  holds  (Appendix  B) 


(3.17) 


o  o 


Substituting  Eq.(3.14)  and  Eq.(3.17)  into  Eq.(3.7),  taking  the  limit  Xp  -►  F 
and  applying  the  boundary  conditions  on  the  cracks,  Eq.(2.3),  the  following 
discretized  BIE's  are  obtained 


in, 

ffl2(Xp,t) 


J  t 
2  ; 
j-l  o 


iri/ 

^22  P  * 


HI21AU1 


j+l  S.  . 

-;j  H‘2iAu  dx^ 

"j  "j 


+  p  /  H|2j^An3^<iXj^ldr 


J  t 
S  / 
j-l  o 


HI22AU2 


J+l 


^  ^222^’^2 ,  l*^^! 


J  j 


s .  - 

+  p  J  ^  H^22^'^2‘^ 


’j 


ijdr  . 


(3.18) 


(3.19) 


The  integrals  in  (3.18)  and  (3.19)  are  to  be  understood  in  the  sense  of 
Cauchy  principal  values.  It  should  be  noted  here  that  Eq.(3.18)  and 
Eq.(3.19)  are  two  decoupled  BIE's  for  the  unknown  crack  opening 
displacements  du^  and  dU2  (as  well  as  their  derivatives) .  The  singular 
terms  of  Eq.(3.18)  and  Eq.(3.19)  at  Xp  -  x  can  be  integrated  analytically  or 
numerically  without  difficulties  as  in  the  usual  "displacement"  BIE 
formulations.  In  the  next  section,  we  will  present  a  time-stepping  scheme 
for  solving  Eqs.(3.18)  and  (3.19). 


To  solve  the  BIE's  (3.18)  and  (3.19),  discretization  of  time  t  is 


necessary.  Here  we  have  used  equal  time  increments  At,  where 
t^  -  nAt  (n  -  1,2,**»N)  denotes  the  time  after  the  n-th  time-step.  The 
unknown  crack  opening  displacements  Au^(x,r)  in  (3.18)  and  (3.19)  are 
approximated  by  the  following  interpolation  functions 


-22  /ij(x)f;"(r)(Au^)j 


(4.1) 


where  (x)  and  »j^(r)  have  the  properties 


(4.2) 


In  Eq.(4.2),  x^  defines  the  i-th  nodal  point,  and 


(Au  )“  -  Au  (x^,r  ) 
a  1  a  —  m 


(4.3) 


represents  the  crack  opening  displacements  at  node  i  and  at  time  mAt. 

In  our  analysis  the  function  (x)  has  been  taken  to  be  unity  over  each 
element  except  for  elements  near  crack  tips.  For  these  elements  a  special 
shape  fvinction 


Mj(x)  -  (a+x^)**  (4,4) 

has  been  used  to  describe  the  proper  behavior  of  Au^  at  the  crack  tips  x,  - 
±  a.  Higher  order  shape  functions  for  »»^(r)  are  desirable  since  Eqs.(3.18) 


Ijl 


and  (3.19)  contain  not  only  the  functions  but  also  their  derivatives 

In  this  paper  the  piecewise  linear  shape  function 


[r)  -  ■ 


1-  ^ ^  '  |r-nAt|  <  At  , 
0  ,  otherwise  , 


has  been  employed. 

For  each  time-step  Eqs.(3.18)  and  (3.19)  can  be  rewritten  as 


in,  ^  . 

‘^la^^Sp-V 


N  J  r 

-  s  s  : 

n-l  j-1  L 


{  l5;5i<5p'5>  if  Mj(x)dx^  + 


p  X  ll^j^(Xp:x)/ij(x)dXj^j(Au^)"  , 


In,  ^  . 

‘^22^5p’ V 


n-l  j-1  *- 
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D^*^(r,q)  Is  defined  as 
mn 


.  f  [(m-n-q)2-r* 

'jo. 


/(c,At)2]P/^  ,  if  (m-n-q)  >  r/(c-At) , 


otherwise . 


(4.13) 


By  choosing  J  collocation  points  on  T  and  N  points  for  t,  and  requiring 

that  Eqs.(4.6)  and  (4.7)  are  satisfied  at  each  discrete  point 

Xp  (i  •  l,2,**oJ),  we  obtain  a  system  of  linear  algebraic  equations  which 

must  be  solved  at  each  time  t  -  mAt  (m  -  1,2,»*»N): 

m 


N  J 


f^-  S  S  A^(Au  )" 
^  n-1  J-1  ’■J  ^  J 


N  J 

gj  "  S  Z  B. ,(AUn).  ^ 

^  n-1  j-1  ^  J 


(4.14) 


(4.15) 


where 


^  “  ‘^12^-P’^m^ 


m  in,  i  _  X 

gi  -  a22(Xp.t^) 


(4.16) 


(4.17) 


-  L^i(Xp:x)/<^(x)  (x)dx, 

^  -^’s.  s.  1-^ 


+  P  /  M^j^(Xp;x)#*.  (x)dXj^  , 

e  J 


(4.18) 


^  ^222^^P’5^i^  Mj(x)dx^ 

s .  s ,  1 


+  P  /  M^2(5p:*)  M.(x)dx^  . 

c  » 


Here  we  will  discuss  only  in  some  more  detail,  since  the  main 
structure  of  and  is  the  same.  At  first  sight  it  seems  that 
discrete  kernels  A^  have  to  be  calculated.  However,  this  number  can  be 
reduced  to  NJ^ ,  if  we  use  the  causality  properties 


‘"a^7^5p-t„:x,t^)  -  0  ,  ifn>m  , 


and  the  following  time  translation  properties  of  the  Green's  functions 
(Appendix  A) : 


-  ‘'5/97^5p' V  ^  ' 


Thus ,  the  matrix  A^  has  the  following  special  form 


where  q  is  defined  by 


q  -  m-n+1 


(4.25) 


The  form  of  b1^  is  similar  to  A^,  and  it  will  not  be  given  here  for 
tj  r  j 

brevity. 

By  considering  this  special  structure  of  and  B™”  we  obtain  finally 
the  following  time*stepping  scheme: 


'(Au^)[  -  fj  , 


(Au,)^  -  i  (Ai  -  V  a52-P'"^^Au,)P] 

^  i  j.l  IJ  J  p_i  J-  It 


(4.26a,b) 


[(Au  -  S  (Bi  )-^  gj 

j-1 


l(AuO?  -  2  fg?  -  V 


(4.27a.b) 
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in  which  (A|j)  and  (®^j)  denote  the  elements  of  the  inverse  matrix  of  A* 


and  at  the  first  time  step.  Also,  q  •  2,3, •••N. 


At  each  time  step  only  one  blockmatrix  A^^  (and  has  to  be 


evaluated.  For  the  first  time  step,  the  inverse  matrix  (A')*^  and 


must  also  be  determined.  Spatial  integrations  in  Eq.(4.18)  and  Eq.(4.19) 
have  been  performed  analytically  for  the  constant  shape  function  (x)  - 
1),  and  numerically  for  the  "crack-tip"  shape  function  (see  Eq.(4.4))  by 
using  an  8-points  Gauss -Jacobian  formula.  In  the  latter  case,  the  singular 
terms  have  been  integrated  analytically  and  numerically. 

5.  Dynamic  Stress  Intensity  Factors 
Once  the  crack  opening  displacements ,  Au^ ,  have  been  calculated  from 
the  time-stepping  scheme  as  described  in  the  last  section,  stress  intensity 
factors  can  be  calculated  by  using  the  following  well-known  relations  [30) 


riqct)  I  _ 


iJ2^ 


4(1-1/) 


fim 


Xj^-*±a  (a+x 


_i _ rhU2(x^,t)' 

lAu^(x^,t) 


(5.1) 


where  "+"  indicates  the  tip  at  x^  -  a  and  "-"  indicates  the  tip  at  x^  -  -a, 
while  V  denotes  Poisson's  ratio. 

In  the  numerical  calculations,  the  incident  wave  was  taken  as  either  a 
plane  longitudinal  wave  of  the  form 


u^”  -  -  (x^+a)sin^-X2COS^] ‘HfCj^t  -  (Xj^+a)sin^-X2Cosd]  (5.2) 


or  a  plane  transverse  wave  of  the  form 


'in*i*!wS 


-  “t 


(Xj^+a)sin^-X2Cos^3 ‘Hfc^t  -  (Xj^+a)sin^-X2Cos^] ,  (5.3) 


where  and  are  the  displacement  amplitudes,  ^  is  the  angle  of 
incidence,  and  H(*)  is  the  Heaviside  function.  The  corresponding  stress 
components  are 


-  Uj^(A+2>x  -  2MCOs2^)H[Cj^t-(Xj^+a)sin^  -  X2Cos^]  ,  (5.4) 


“  UL/<sin(2^)H[Cj^t-(Xj^+a)sin^  -  X2Cos^]  ,  (5.5) 


<^22  “  ■  2p»sin2^)H[Cj^t-(Xj^+a)sin^  -  X2Cos^]  ,  (5.6) 


for  an  incident  longitudinal  wave,  and 


‘^ll  “  *U^sin(2^)H[c,j,t-(x^+a)sin^  -  X2Cos<^]  ,  (5.7) 


<^12  “  ■  U.j./iCos(2^)H(c.j,t-(Xj^+a)sin^  -  X2Cos^] ,  (5.8) 


<^22  "  U^Msin(2^)H[c.j,t-(Xj^+a)sin^  -  X2Cos^]  (5.9) 


for  an  incident  transverse  wave. 

All  calculations  have  been  carried  out  for  a  Poisson's  ratio  i/  -  1/4. 
The  geometrical  configuration  is  shown  in  Fig.  1.  The  principal  (macro) 
crack  has  been  discretized  into  50  elements  of  equal  length,  and  a 
proportional  number  of  elements  have  been  used  for  the  micro-crack. 
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For  b/a  -  0,  the  configuration  reduces  to  a  single  crack  of  length  2a. 
This  case  was  used  to  check  numerical  results  obtained  by  our  method. 
Comparisons  with  Thau  and  Lu's  results  [2]  are  given  in  Figures  2a  and  2b. 
Figure  2a  presents  results  for  normal  incidence  of  a  longitudinal  (L)  wave, 
while  Fig.  2b  is  for  a  normally  incident  transverse  (TV)  wave.  All  results 
are  normalized  by  the  corresponding  static  values.  The  agreements  are  very 
good,  especially  in  the  overall  variation  of  the  stress  intensity  factors 
with  time  and  in  their  peak  values.  Calculations  for  inclined  incidence 

0)  of  L  or  TV  waves  have  also  been  carried  out,  and  the  agreements  with 
Thau  and  Lu's  results  are  again  very  good.  The  time  Increment  was  selected 
as  c.j^t  -  0.08a.  The  influence  of  At  on  the  stability  of  the  time-stepping 
scheme  has  been  studied  numerically,  and  it  was  found  that  too  small  a  value 
of  At  may  cause  instabilities  for  the  results  at  large  time.  The  same 
conclusions  have  been  drawn  by  Nishimura  et  al.[13].  The  time  increment 
chosen  here  always  yielded  good  results,  and  hence,  this  value  has  also  been 
used  for  the  macro-microcrack  interaction  problems. 

For  a  fixed  half-length  of  the  microcrack,  b/a  —  0.1,  and  for  normal 
incidence  of  a  longitudinal  wave  -  0°),  the  dynamic  stress  intensity 
factors  are  shown  in  Figures  3a, b,  versus  the  dimensionless  time  c^^t/a,  for 
various  values  of  the  crack-tip  distance  d/a.  All  results  have  been 
normalized  by  the  static  stress  intensity  factors  of  a  single  crack  under 
the  corresponding  static  load.  Due  to  symmetry  with  respect  to  X2  -  0,  the 

Mode-II  stress  intensity  factors  are  identically  zero,  -  0.  Figure  3a 
shows  that,  as  expected,  the  presence  of  a  microcrack  does  not  influence  the 
left  tip  of  the  macrocrack  at  small  time.  However,  after  c^^t/a  =«  3  the 
difference  between  both  cases  becomes  somewhat  distinct.  It  is,  however. 


evident  that  the  crack- tip  away  from  the  microcrack  is  not  significantly 
affected  by  the  presence  of  a  colllnear  microcrack.  The  presence  of  a 
microcrack,  does,  however,  give  rise  to  a  substantial  increase  of  the  stress 
intensity  factor  at  the  tip  adjoining  the  microcrack  (Fig.  3b).  The 

K^-factor  increases  with  decreasing  crack-tip  distance  d/a,  and  considerable 

amplifications  in  occur  for  very  small  values  of  d/a.  At  both  tips  of 
the  main  crack,  the  maximum  dynamic  stress  intensity  factors  exceed  the 

—  I 

corresponding  static  values,  which  are  reached  at  large  time  t.  Both  and 

K^,  assume  the  values  for  a  single  solitary  crack  of  half-length  a  as 
d/a  -*  ®. 

Figures  4a, b  and  Figures  5a, b  show  the  normalized  dynamic  stress 
intensity  factors  when  a  longitudinal  wave  is  indcident  under  an  angle 
— “f-  — <4* 

<i>  -  30® .  Both  the  and  factors  are  present  in  this  case .  The  time 

history  of  Kj  (Fig.  4a)  is  very  similar  to  that  for  normal  incidence 

(Fig.  3a).  The  K^-factor  at  the  tip  adjoining  the  microcrack,  which  is  zero 
before  the  incident  wave  front  arrives  at  this  tip  (c^t/a  <  1) ,  is  shown  in 
Fig.  4b.  The  smaller  the  distance  between  microcrack  and  main  crack,  the 
larger  is  the  stress  intensity  at  the  crack-tip  adjoining  the  microcrack. 

It  is  also  interesting  that  the  maximums  of  the  normalized  values  are  the 

same  as  those  for  normal  Incidence.  Results  for  are  given  in  Figures 
5a, b.  Their  variations  with  the  dimensionless  time  c^t/a,  and  their 
dependence  on  the  crack- tip  distance  d/a,  are  similar  to  those  of  the 
— + 

K^- factors.  The  overshoots  of  the  Mode- II  d}mamlc  stress  intensity  factors 
are,  however,  less  than  those  of  the  Mode- I  cases. 


It  is  evident  that  the  dynamic  stress  Intensity  factors  depend  strongly 
on  both  the  location  and  the  size  of  the  microcrack.  For  a  fixed  crack- tip 


distance,  d/a  —  0.05,  and  for  normal  incidence  of  a  longitudinal  wave 

—  ^ 

-  0°),  the  dependence  of  the  K^- factors  on  the  dimensionless  half-length 
of  the  microcrack,  b/a,  is  shown  in  Figures  6a, b.  As  before,  at  small  time, 
the  crack- tip  away  from  the  microcrack  behaves  as  the  tip  of  a  semi- infinite 

crack,  and  Kj  is  the  same  for  all  five  cases  (proportional  to  .  The 

influence  of  the  microcrack  (as  well  as  the  crack- tip  adjoining  the 

microcrack)  on  becomes  important  when  diffracted  waves  arrive  at  the  left 

tip  of  the  main  crack.  The  peak  K^-factor  increases  with  increasing  b/a, 

but  it  is  shifted  to  somewhat  larger  time  Cj^t/a.  The  K’ -factor  for  larger 
b/a  can  be  slightly  smaller  than  for  a  shorter  microcrack  (smaller  b/a), 
which  is  in  contrast  to  the  static  case  (see  Yokobori  et  al.[25]).  Figure 

6b  shows  the  variation  of  Kj  with  che  microcrack  size  b/a.  As  expected,  a 
larger  microcrack  gives  rise  to  a  larger  amplification  of  the  stress 
intensity  factors.  For  b/a  -  0  the  configuration  reduces  to  a  single  crack, 

* 

and  in  this  case  we  have  • 
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The  Green's  Function 

The  two-dimensional  elastodynamic  Green's  function  is  given  by  (see 
[281.  [29]) 
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where 


r  -  |x-Xp|  ,  -  yc|(t-r)2.r2,  ^  -  L.T  ,  (A. 2) 

and  H(»)  is  the  Heaviside  step  function.  The  function  u“^(Xp, t;x, r)  denotes 
the  displacement  in  the  a-direction  observed  at  position  x  and  at  time  t, 
due  to  a  unit  force  in  the  7  direction,  applied  at  position  Xp  and  at  time 
r.  The  corresponding  stress  components  follow  from  (A.l)  and  Hooke's  law  as 
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K  -  Cj^/C^  . 

The  function  u^^(Xp,C;x,r)  has  the  following  properties 

1)  causality 

'"a7^5p’*^  =  -’’‘^"  °  ^  l5'5pl  ’ 

2)  time  translation 
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translation  properties . 


(Xp,t;x,T)  .  (A. 9) 

also  possess  these  causality  and  time 


To  prove  the  identity,  Eq.(3.17),  it  is  sufficient  to  consider  only  the 


first  term  of  since  the  three  terms  of  H*.,.  have  similar  forms.  We 
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Partial  integration  of  the  right  hand  side  of  (B.l)  yields 
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The  first  two  terms  in  (B.2)  vanish  due  to  the  initial  conditions  on  Au^  and 


G  *0 

AUj5  ,  and  due  to  the  causality  properties  of  u^^  and  u^^  .  Thus ,  we  obtain 


S  -  X  u®  Ad^dr 

o'  o' 


(B.3) 


Equation  (3.17)  can  be  derived  directly  from  (3.16)  and  (B.3). 
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Fig.  2:  Normalized  dynamic  stress  intensity  factors  for  a  single  crack, 

and  for  normal  incidence  of  a  plane  wave, -  Thau  and  Lu, 

o  o  o  present  work;  (a)  Kj  for  L-wave  incidence;  (b)  for  TV 
wave  irrcidence. 
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Normalized  stress  intensity  factors  for  various  d/a,  and  for  normal 
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incidence  of  a  plane  L-wave,  b/a  =  0.1;  (a)  K,;  (b)  K,. 


f*  ••*  •»*»  1 


Unclassified 


security  Classification  of  this  pace  (Wh«n  nat*  Fntorad) 


REPORT  DOCUMENTATION  PAGE  BEFolE'co^^^^^NrPORM 


I.  REPORT  number  2.  GOVT  ACCESSION  NO.  3  RECIPIEN  T'S  C  AT  alOG  NUMBE  R 

NU-SML-88-1 


«.  TITLE  Subtlll*) 


Time-Domain  Boundary  Element  Analysis  of 
D)naamic  Near-Tip  Fields  for  Impact-Loaded 
Collinear  Cracks 


7.  AUTHORCa; 


Ch.  Zhang  and  J.  D.  Achenbach 


5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Interim 


6  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERfaj 


N00014-85-0401 


9  PERFORMING  OROANI  ZATION  NAME  AND  ADDRESS 


10.  PROGRAM  element.  PROJECT  TASK 
AREA  A  WORK  UNIT  NUMBERS 


Northwestern  University,  Evanston,  IL  60208 


II.  controlling  OFFICE  name  ANO  AOORESS  12.  REPORT  DATE 

Office  of  Naval  Research  February  1988 

Structural  Mechanics  Department  number  of  pages 

Department  of  the  Navy,  Arlington,  VA  22217  32 


14.  monitoring  agency  name  a  AODRESSf/f  dllittmt  (tom  Comroitint  OfCicaJ  15.  SECURITY  CLASS,  (at  Ihia  report) 

Unclassified 


'5»,  DECL  ASSI  FICATION/ DOWNGRADING 
SCHEDULE 


16.  distribution  statement  (ot  Ihle  (iapart) 


Approved  for  public  release;  distribution  unlimited 


17.  OISTniBUTION  STATEMENT  (a/  ttf  abatract  antarad  in  Biock  20,  if  diffarent  from  Raport) 


19.  KEY  WORDS  (Continua  on  ravaraa  stda  if  nacaaaary  and  fdantify  by  block  numbar) 

elastodynamic  stress  intensity  factors 

macrocrack 

microcrack 

time-domain  boundary  element  method 


20.  abstract  (Continua  on  ravaeaa  aida  II  nacaaaary  and  ldant*fy  by  block  numbar) 

A  time-  domain  boundary  integral  equation  method  has  been  developed  to 
calculate  elastod)uiamic  fields  generated  by  the  incidence  of  stress  (or 
displacement)  pulses  on  single  cracks  and  systems  of  two  collinear  cracks.  The 
system  of  boundary  integral  equations  has  been  cast  in  a  form  which  is  amenable 
to  solution  by  the  boundary  element  method  in  conjunction  with  a  time-stepping 
technique.  Particular  attention  has  been  devoted  to  dynamic  overshoots  of  the 
stress  intensity  factors.  Elastodynamic  stress  intensity  factors  for  pulse 
incidence  on  a  single  crack  have  been  computed  as  functions  of  time,  and  thev 


1  JAN  71  1473  edition  of  I  NOV  65  IS  OBSOLETE 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  thiS  PAGE  (When  Dai*  Enterod) 


have  been  compared  with  results  of  other  authors.  For  collinear  macrocrack-micro 
crack  configurations  the  stress  intensity  factors  at  both  tips  of  the  macrocrack 
have  been  computed  as  functions  of  time  for  various  values  of  the  crack  spacing 
and  the  relative  size  of  the  microcrack. 


])Ar£L 

piL  b2> 

Dr/c- 


